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Abstract 

Passive random walker dynamics is introduced on a growing surface. The walker is designed 
to drift upward or downward and then follow specific topological features, such as hill tops or 
valley bottoms, of the fluctuating surface. The passive random walker can thus be used to directly 
explore scaling properties of otherwise somewhat hidden topological features. For example, the 
walker allows us to directly measure the dynamical exponent of the underlying growth dynamics. 
We use the Kardar-Parisi-Zhang(KPZ) type surface growth as an example. The word lines of a set 
of merging passive walkers show nontrivial coalescence behaviors and display the river-like network 
structures of surface ridges in space-time. In other dynamics, like Edwards- Wilkinson growth, this 
does not happen. The passive random walkers in KPZ-type surface growth are closely related to 
the shock waves in the noiseless Burgers equation. We also briefly discuss their relations to the 
passive scalar dynamics in turbulence. 

PACS numbers: 05.40.-a, 02.50.Ey, 64.60. Cn, 89.75.Da 
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I. INTRODUCTION 



The research on fluctuations of surfaces during growth has been one of the major areas 
of study in nonequilibrium statistical mechanics in recent decades. Most research is focused 
on identifying universality classes and the scaling behavior of surface morphology, e.g., how 
the width of surfaces scales with the size of surfaces [| |], ||, |[] . Recently, the coupling 
between surface growth and other degrees of freedom has been considered in the context of 
ordering phenomena on growing and fluctuating surfaces. For example, the field coupled to 
the surface can be different species of deposited particles^], ^, |8| or surface reconstruction 
order parameters P]. Unlike equilibrium surfaces, where the coupling between the surface 



height degrees of freedom and the additional field is irrelevant at large length scales [10 



we found that the coupling between the surface and the reconstruction orders can not be 
ignored 0. The surface and the additional field are usually coupled topologically. For 
example, the domain walls of order parameters become trapped at the hilltops or in the 
valleys on surfaces in 1+1 D system^, [7], || or on the ridge lines on surfaces in a 2+1 D 
system 0. 

The fluctuations of the additional degrees of freedom (expressed by the dynamics of the 
"domain walls") thus become slaved to the fluctuations of the surface. From the surface 
perspective, the additional field is usually irrelevant and passive. Therefore, part of the 
dynamics of the surface and its topological features can be evaluated from the fluctuations 
of such domain walls that are pinned or trapped to them. We can put probes on the surface 
to follow these features dynamically. In this paper, passive random walkers(PRWs) are 
designed to follow hilltops or valley bottoms on the surface. 

One of the first applications of the PRWs is to provide a direct way to measure the 
dynamical exponent of the surface fluctuations. We will also present the phenomena of 
coalescence of passive random walkers on Kardar-Parisi-Zhang(KPZ) Jl|, type growth. The 
coalescence of passive random walkers uncovers a river-like network structure of the surface 
in space-time. 

Our passive random walkers are similar to the passive scalars in turbulence. This is not 
surprising because the KPZ equation is equivalent to Burgers equation, which describes a 
fluid |1| with different random stirring forces from those in related turbulence studies. It is 
instructive to study the dynamics of the PRWs from both perspectives. 
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This paper is organized as follows. The passive random walkers model is defined in Sec. II. 
The application to directly measuring the dynamical exponent is presented in Sec.fT|. Then 
we discuss the coalescence of the PRWs in Sec.[TV]. The sign of the coupling between the 
PRWs and the surface is important for the coalescence phenomena as discussed in Sec.|VT|. 
In Sec. [VTTL we show the river-like network in the space-time structure of KPZ-type surfaces 



and the relation between the passive random walkers and the passive scalars in turbulence 
is briefly discussed. Other possible applications of the PRW model are pointed out, together 
with a summary of the results in Sec. [VIII| . 

II. PASSIVE RANDOM WALKERS ON A GROWING SURFACE 

We couple a PRW to the growth dynamics of a fluctuating surface. The movement of the 
PRW is determined by the local slope of the growing surface. The well-known Kim-Kosterlitz 
(KK) model0 for surface growth is used as an example. In the KK model, a surface is 
specified by integer height variables on a square lattice. A single particle is deposited at a 
randomly chosen site i if the restricted solid-on-solid condition = \hj — hi\ < 1 for 

all nearest neighbor pairs (ij)) remains satisfied, otherwise, the deposition of the particle 
is rejected. The surface grows and the stationary state is rough. The position of the PRW 
is updated after N^/v g Monte Carlo steps (with v g the growth velocity and N s the total 
number of site), i.e., when on average one layer of surface material is added, according to 
the following rules. Let i be the position of the walker. If there is any neighbor site j for 
which hj > hi, then we move the walker to site j. If there is more than one site higher than 
site i, then we let the PRW move to either one with equal probability. We apply periodic 
boundary conditions for both the surface and the PRW motion. 

It is known that the continuum limit of the KK model is governed by the so-called KPZ 
equation 0, 
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dh d 2 h , x (dh\ 
with uncorrelated noise, 

(r) s (x,t)r) s (x',t')) = D6(x,x')8(t,t'). (2) 

In the same limit, the equation of motion of the PRW is 

du d 



where u(t) is the coordinate of the PRW at time t, and k is the coupling strength between 
the PRW and the surface. In our lattice model, k is of the same order as A. The equations 
of motion, Eqs.([I])-(|3]), imply that the PRW moves upwards for k > 0. The PRW becomes 
trapped on a local maximum of the surface profile as seen in Fig. (|Tj) for positive n. Once 
the PRW is trapped on a hilltop, it follows the motion of the hilltop, which is subject to the 
surface fluctuations. When k is negative, the PRW moves downwards instead of upwards, 
and the PRW becomes trapped in a valley bottom instead of a hilltop. The relative sign 
between k and A affects the dynamics of the PRWs. We will focus on positive k for now, 
and discuss negative n in Sec.|V]]. 



III. DYNAMICAL EXPONENT OF PASSIVE RANDOM WALKERS 

One important observation of the dynamics of an upwards moving PRW on KPZ-type 
surfaces is that the PRW not only reaches a local maximum in a relatively short time, but 
also that in finite-size systems it ultimately ends up trapped at the global maximum. We 
will discuss the detailed mechanism of this phenomenon later in Sec.^. Once the PRW is 
trapped on the global maximum, it performs a correlated random walk following the global 
surface fluctuations. The surface fluctuates critically and is self-affine. Namely, the profile 
of the surface is invariant under the following rescaling, 

x -> b^x, h -> b' a h, t -> b- 1/z H, (4) 

where b is a scaling factor, a is the roughness exponent, and z s is the dynamical exponent 
of the surface. The PRW is slaved to the fluctuations of the interface, therefore we expect 
similar dynamical scaling for the world line of the PRW. 

The displacement of the PRW, Au(£) = u(£) — u(0), should be defined carefully to avoid 
confusion in finite-size systems with periodic boundary conditions, where the PRW moves 
on a ring in ID and on a torus in 2D. In such manifolds, the winding number of the PRW 
should be taken into account. The component of the displacement, Au(t), in direction ej is 
defined as Ait§.(t) = nf. — r%, with the actual number of right moves, raj, and the actual 
number of left moves, in direction e^. 

In our numerical simulations, we let the system evolve until both the surface and the 
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FIG. 1: Passive walkers trapped on local hilltops and global maximum in 1+1 D simulations. The 
dashed lines are the world lines of two PRWs and the solid lines are the surface profiles during 
growth. The passive walkers follow the local maximums initially and ultimately merge into each 
other. The walkers reach the global maximum at a later time. 
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PRW motion reach stationary states. Then the value of the displacement Au(t) is measured 
as Au(t) = (nti ~ n h) ) - To accelerate the simulation, we also adopt a rejection- 
free algorithm. Details of this algorithm can be found in ref. 0. In this algorithm, time is 
counted in terms of the numbers of layers of particles deposited instead of the conventional 
Monte Carlo time unit. We established earlier that the unit of time in this rejection free 
algorithm is linearly proportional to the Monte Carlo time unit (see also ||). 

If the PRW indeed is slaved to follow the global surface fluctuations, which are invari- 
ant statistically under the transformation eq.([|), the average distance of the displacement, 
Au(t), must obey the scaling form, 

<Au(f)> ~ L X Q (±) , (5) 



where Q{r) is a universal scaling function. When t <C L Zs , the hopping events of the PRW 
are correlated due to the surface fluctuations, and Au(t) ~ t Zw , where z w is the dynamical 
exponent for the PRW. At time scales t 3> L Zs , the surface fluctuations are limited by 
the finite size of the lattice and the PRW becomes like a free uncorrelated random walker. 
Thus, the scaling function Q{t) has the asymptotic forms, G(t) — > r 1 / 2 ™ for r ■€ 1 and 
Q(t) — > r 1 / 2 for r ^> 1. The value of the exponent % in Eq.(|5|) follows from the fact that 
u(t) must be independent of L in the thermal dynamics limit. Consequently, \i z w, and z s 
are not independent and satisfy the relation \ = jj*-. 

The dynamical exponent of the walker is not necessarily the same as the dynamical 
exponent of the surface. There is no a priori reason for % = 1, i.e., for z w = z s . The 
numerical results (Fig. |2|) indicate that z w = z s , but why? \ = 1 reflects that there is no 
length scale other than the system size L involved in the dynamics of the PRW. This is in 
contrast to free uncorrelated random walks, in which case another length scale proportional 
to if: 1 / 2 will emerge. The absence of a new independent length scale is consistent with our 
picture that the PRW is typically trapped on a global maximum, hence the fluctuations of 
the PRW follow the dynamical scaling of the surface. 

The main numerical results are shown in Fig.(^)-([|). They are obtained by averaging over 
about 10 5 samples in 1+1D and 4 x 10 5 samples in 2+1D. The displacements of the PRW 
are determined up to t = 256 in 1+1D and t = 128 in 2+1D. In 1+1D, (Ait) is in the order 
of 10 2 lattice units at t = 256. In 2+1D, (An) ~ 10 at t = 128. Although (Au) is relatively 
small in 2+1D, the scaling behaviors are already clear. Small displacements reflect merely 
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a small prefactor in the scaling form eq.(j5|). This is analogous to a small diffusion constant 
(which slows the diffusive process, but does not affect the scaling properties even at small 
spatial scales). Moreover, the asymptotic region of interest is at small r and at small Au/L, 
so the scaling behaviors that we focus on are in the small Am interval. By increasing system 
sizes, we can extend the scaling region of the scaling variable, Au/L, about two decades. 

We check the scaling form eq.(|^) by collapsing the data (see Fig. (0)) with z s = 1.5 in 
1+1 D and z s = 1.6 in 2+1 D. The prediction, \ — 1> gives excellent data collapse, except 
for the nonuniversal part around Am < 10 in 1+1D and Am ^ 4 in 2+1D. This nonuniversal 
behavior is likely caused by the discrete lattice spacing. A more detailed fit of the dynamical 
exponent z w must avoid this small Au region while keeping Au/L small. In Fig.(|3|b)(P3), we 
fit z w by omitting the data from Au < 10 in 1+1D and Au < 4 in 2+1D. It might appear 
that by ignoring data from these regions, a large fraction of information in the data is lost. 
In fact, only less than a quarter of the data points fall in these regions. The crossover from 
uncorrelated random walker behavior (z = 2, at large Au/L) to correlated behavior (z ^ 2, 
at small Au/L) is visible in Fig.(|3|b)(|4]b), where the exponent decreases from 2 to z s as the 
system size increases. We conclude that z w = 1.50(1) in the 1+1D and z w = 1.60(1) in the 
2+1D KK model. 

Our measurement of z w provides one of the very few independent and direct measure- 
ments of the dynamical exponent z s on stationary- state KPZ-type surfaces. Conventional 
approaches determine the dynamical exponent indirectly, through the measurement of the 
scaling behaviors of the width of surfaces in the transient states, starting with a flat initial 
condition. The surface width defined as W = ((h—h) 2 ) 1 ^ 2 scales as W ~ t 13 for t <C L Zs . Then 
the dynamical exponent is found by using the scaling relation z s = a/f3. The other method 
is to measure the correlation function g(x, t) = ((h(xo + Ax, t + At) — h(xo, to)) 2 ) ~ (At) 2/3 
for At < (Ax) 2s at stationary states. Both methods obtain the dynamical exponent indi- 
rectly through the growth exponent (3. By tracing the path of the PRW, we are able to 
provide an independent and direct measurement of the dynamical exponent without knowl- 
edge about the global width of the surface. 

Compared to the conventional methods of measuring the dynamical exponent, we only 
need to acquire information of the local slopes around the PRW. Thus, it is instructive to 
understand the connection between the slope-slope correlations along the path of the PRW 
and the global roughening dynamics of the surface. 



7 




FIG. 2: Collapses of t versus Au. We use z s = 1.5 for 1+1D, z s = 1.6 for 2+1D and % = 1 for both 
cases. 
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FIG. 3: (a)The log-log plot of t versus An in the 1+1D KK model. For reference, we plot the solid 
line with a slope of 1.5. (b)Fitting of z w for different system sizes in 1+1D. It shows the crossover 
from free uncorrelated random walker behavior (z = 2) to a KPZ-type PRW behavior (z = 1.5). 
The dotted line is a convenience to guide the eyes. 
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FIG. 4: (a)The log-log plot of t versus An in the 2+1D KK model. For reference, we plot the solid 
line with a slope of 1.6. (b)Fitting of z w for different system sizes in 2+1D. It shows the crossover 
from free uncorrected random walker behavior [z = 2) to a KPZ-type PRW behavior (z = 1.6). 
The dotted line is a convenience to guide the eyes. 
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The displacement of the PRW and the correlation of the slopes are connected by the 
fluctuation-dissipation theory. The scaling of the dispersion in the PRW displacement can 
be evaluated from the slope-slope correlation function along the world line of the PRWflTT 
i.e., 



with $(V) 



((u(t) - u(0)) 2 ) ~ t f dr $(r) ~ t 1/z ™, (6) 
Jo 



This correlation function <&(r) is different from the conventional slope-slope correlations, 
$(r, t) = (d x h(r, t)d x h(0, 0)), where the correlations are calculated at specific fixed r. In- 
stead, we need to evaluate the correlation $(r, t) along the correlated path, the world line 
of the PRW, r = u(t). 

The most general scaling form for $(r, t) is, 



$(r,t) = 6 2 "l> (W,r 2s t) , (8) 
where rj is an exponent yet to be determined. This scaling form is equivalent to 

$(r,t) ~t 2ri/z *F (r-H) , (9) 



with F(y) a scaling function. Since a path of the PRW is determined by the surface fluc- 
tuations, averaging over different paths can be approximated by averaging over different 
realizations of the surface (like different Monte Carlo runs). $(t) is approximately equal to 
the value of $(r, t) at r = Au(t) ~ t 1 ^™. Thus, $(t) scales as 

$( T ) ^ ^(r 1 / 2 ", r) ~ T 2ri/Zs F (r 1 ^) . (10) 

If X > 1, asymptotically, the scaling function behaves as F (r 1_x ) ~ F(0) + F'(0)r 1_x , 
which approaches a constant, for r — > oo. Hence the correlation function scales as $(r) ~ 
T 2r] / Zs . Together with Eq. (f|), we obtain 2 + 277/2,, = 2/^ w . This relation allows us to 
determine 7], which is the scaling dimension of the slope operator, from the measurement 
of z w . Furthermore, if % = 1, then rj = a — 1, provided that the KPZ scaling relation 
a + z s = 2 holds. This result, rj — a — 1, is consistent with Lassig's operator production 
expansion scheme Q and our previous work0 for KPZ type surfaces. Namely, all slope- 
slope correlations can be obtained by the naive power counting and the slope operator, ||, 
scales as x a ~ l . 
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Our numerical results suggest that x is indeed equal to 1 for upwards moving PRWs. 
In 1+1D, renormalization group (RG) calculations gave so our 

results imply the value of z is fixed as z w = z s = 3/2 in 1+1 D. In 2+1 D KPZ-type 
surface growth, RG calculations indicated that the fixed point is at strong coupling such 
that analytical results can not be obtained perturbatively|I|]. To author's knowledge, no 
analytic result for the scaling from of $(r, t) has been obtained for 2+1 D. Our numerical 
results support r\ = a — 1 for KPZ-type surfaces. 

IV. COALESCENCE OF PASSIVE RANDOM WALKERS 

Why does an upwards moving PRW typically end up on the global maximum? To further 
understand the mechanism of this phenomenon, we study the dynamics of multiple PRWs 
on KPZ-type surfaces. If all upward moving PRWs eventually reach the global maximum 
on the finite-sized KPZ-type surface, then all of them will converge to the same hilltop even 
though they start at different positions. This is true only if the PRWs have the property of 
coalescence, namely two PRWs close to each other merge into each other and move together 
afterward during the process of moving up to the global maximum. 

Although, in the definition of the dynamics of the PRWs, the PRWs always move upward, 
this simple coupling between the PRW and the surface only guarantees that the PRW reaches 
a local maximum, not the global one. However, we find that in KPZ dynamics with negative 
A the upwards moving PRWs coalesce. Two PRWs separated over a distance r at time t = 
merge into each other after a typical time period t ~ r Zw as shown in Fig.(^J) ( see also (|8]a) 
). On the contrary, downwards moving PRWs do not coalesce in A < KPZ-type dynamics 
and PRWs do not coalesce in Edwards- Wilkinson(EW) type growth (the A = point of 
KPZ equation). 

Consider the PRWs sitting on two nearby hilltops that have similar heights. The region 
between the two hilltops is a small and shallow valley. If both the hilltops and the valley are 
relatively higher than the other parts of the surface, the valley is likely to be filled up and 
vanish in a short time. These two hilltops merge and the merging event causes the PRWs 
to coalesce afterward as shown in Fig. ([!]). The coalescence of the PRWs is likely caused by 
merging of two nearby hilltops by filling the valley in between. 

This trivial argument seems to imply that PRWs will coalesce for any type of growing 
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surface as long as the PRWs are trapped on hilltops. This is not true. We found that the 
PRWs on EW type surfaces do not coalesce, although KPZ-type surfaces and the EW-type 
surfaces have the same stationary state and roughness exponent a = 1/2 in 1+1 D. The 
following discussion and the detail analysis in Sec.|V] show how the nonlinear term, A (d x h) 2 , 
in the KPZ equation(eq.([I|)) affects the phenomenon of coalescence. 

The nonlinear term, X(d x h) 2 , in the KPZ equation controls how the growth rate depends 
on local slopes. Consider a plateau on a surface, for positive A, the growth rate on the 
slope parts of the plateau is larger than the flat top. The flat top of the plateau expands 
by lateral growth on the slopes. In general, hilltops become flatter than valleys for positive 
A. Negative A has an opposite effect, the slope parts grows slower than the flat tops hence 
the flat tops become sharper. Shrinking of plateau flat tops is an indication of a negative 
A in the KPZ equation. The merging events of two hilltops is equivalent to shrinking of 
flat top plateaus. For example, at larger length scales, the shallow valley between two 
nearby hilltops and the hilltops themselves as a whole acts like a plateau on the surface 
(Fig. (|5])). Before the coalescence of the hilltops, the size of the flat top of the plateau is 
the distance between these two hilltops. While these two hilltops are moving towards each 
other by surface fluctuations, the size of the flat plateau top effectively shrinks, which is an 
indication of a negative A. In the KK model, A is negative and the origin of negative A is 
due to the restricted solid-on-solid condition |A/ty| < 1, which reduces the growth rate to 
zero at the regions where the density of Ah = l(or Ah = —1) steps is 1. With negative A 
for the KK model, which belongs to the KPZ universality class, both the hilltops and the 
PRWs coalesce. 

V. COALESCENCE OF HILLTOPS IN THE NOISELESS KPZ EQUATION 

We can gain insights about the coalescence of the PRWs on KPZ-type surfaces by inves- 
tigating the solution for the noiseless KPZ equation in the limit v — > 0. In that limit, the 
evolution of surfaces with given initial conditions has an analytical solution, which allows 
us to study the details of the dynamics of the hilltops. 

A detailed derivation of the solution for the noiseless KPZ equation can be found in 
|T3| . We only present a brief review of the derivation. With the Hopf-Cole transformation, 
(f>(x,t) = exp (Xh(x , t) / u) , we can transform the noiseless KPZ equation into a simple linear 
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FIG. 5: A close look at the coalescence of two hilltops driven by the noiseless KPZ equation. The 
curves in this figure have been shifted for clarity. The two hilltops have a common contact point 
at the bottom of the valley between them. 
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diffusion equation for <f>(x,t), 

d<p d 2 (j) 

The solution is a superposition of Gaussian functions weighted by the initial conditions, 

<f>(x,t) = (2tt^)-' Jdy<f>{y, 0)exp {- ^Z^D j , ( i2 ) 

where <j)(y,0) = exp (Xh(y, 0)/V) is the initial condition, determined by the initial height 
profile 0). In the limit z/ — > 0, we can simplify the solution by introducing the velocity, 
v = d x h, as an auxiliary variable, 

v _ u d In $ _ J <fr ggg exp (-1 + A%, 0))) 

& " /^exp(-i(^ + A%,0)' N 



By the method of steepest descent in the limit v — ► 0, this yields 



t) = 2(X ^fellli , (14) 
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where y min is a function of both x and t, and it is determined by taking y where — 
Xh(y, 0) is a minimum for fixed x and fixed t. Namely, 

f(x-y) 2 \ 
ymm{x,t) = argmin I Xh(y,Q)\ . (15) 

The y m i n (x, t) are usually only piecewise continuous functions of x, particularly for random 
initial conditions generated by random walkers. Shock waves of the velocity field are formed 
where y m i n (x,t) is discontinuous. 

Assuming y m i n (x,t) is constant between x and y m i n , we can further integrate v(x,t): 

h[x, t) = — + c h , (16) 

where Ch is the integration constant. The solution for h(x,t) is reduced to the problem of 
finding y min and c h . 

Determining y m i n for A < is equivalent to finding y for which the vertical distance 
is the shortest between the initial height configuration, h(y,0), and the parabola, ^ x ~ x f , 
centered at x, while the parabola is approaching the surface from below. y m i n is the first 
contact point between the two curves h(y, 0) and *&(y; x, t) = + c while c is adjusted 

to c(x, t) = % m i n , 0) — (Fig. (H)). The contact point y m i n is where h(y, 0) — is 

a minimum. Once j/ min is found, h(x,t) follows from eq.flTB]) with Ch = h(y min ,t). Thus, we 
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have h(x,t) = *f?(x;x,t) = c(x,t), which is the top of the parabola. The surface between x 
and y m in(x, t) is also parabolic in shape at time t. It is possible that there exists more than 
one first contact point for certain x and t. If this happens, y m i n (%, t) as a function of x will 
be discontinuous at that specific x. 

This solution has an intuitive graphical interpretation, which is useful for understanding 
the coalescence phenomena. The evolution of the surfaces is generated by a sequence of 
geometry transformations on the surfaces at an earlier time. For negative A, the new surface 
at time t is obtained by scanning the original surface from below with a probe that has a 
parabolic shaped tip, ty(y;x, t) = ^-^r — l~ c i x i ^)(Fig-®) f° r negative A. (For positive A, 
the probe scans the surface from above instead.) The probe is centered at position x, and 
c(x, t) (the vertical position of the top of the tip) is adjusted by moving the probe up and 
down. When the probe approaches the surface from below, it stops when it starts to come 
into contact with the surface. The horizontal position of the first contact point between the 
probe and the surface is y m m(%,t). Since the probe stops moving upward once it comes in 
contact with the surface, we might think of the vertical position of the tip of the probe, 
c(x,t), as the height of the new surface seen by the probe at x at time t. The probe scans 
through the surface and plots a new surface that gives the surface evolved by the noiseless 
KPZ equation at a later time. Because the parabolic shape of the probe, the new surface 
usually consists of a collection of parabola shaped segments(the dashed curves shown in 
Fig- (I))- 

With this geometrical interpretation of the exact solution of the noiseless KPZ equation, 
we are able to visualize how the singular hill tops are formed and to discuss the dynamics of 
coalescence of the hill tops. If the local average curvature over a region such as between A 
and C in Fig. (||) is larger than the curvature of the tip of the probe at time t, then the tip 
cannot reach the surface everywhere between point A and C. When the probe scans through 
this region, at point B in Fig.(^), there are two first contact points for the parabola, and 
2/min(^) is discontinued at this point. At points where the function y m i n (x) is not continuous, 
the path of the tip forms cusps, which is are hilltops of the new surface. The position of 
the hilltop in Fig.|6| is determined by the positions of the two contact points A, C and the 
curvature of the probe tip. 
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FIG. 6: The parabolic shape probe tip scans through point x a ,Xb, and x c from left to right. Curves 
(a),(b),and (c) are parabola which are centered at x a , Xf,, and x c respectively. Curve (a)/(b) only 
contact with the initial random surface at point A/B. Curve (b) contacts with the initial surface 
at both point A and C. While the probe is scanning through x^, the first contact point jumps from 
A to C. Thus, Uminix) is discontinued at point xi, and the path of the tip creates a cusp at point 
C. The dashed curve is the surface generated by the path of the probe tip. 
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We can also write the solution of h(x, t) as a functional transformation T defined as 

%{h{x, t )) = h(x, t + t )= mm (h(y, t ) - ^ ~ V) ^ . (17) 

This transformation has the following property, 

T tl+t2 (h(x,t )) = T t2 T tl (h(x,t Q )). (18) 

The geometrical interpretation of the solution in this form is similar to the Huygens principle 
growth algorithm except that parabola shaped wave fronts instead of circular or spherical 



wave fronts are used. The transformation eq.([PTD can be applied iteratively and is suitable 
for evaluating the evolution of surfaces numerically. 

We evaluate the evolution of a random initial surface with Eq.([T7|). The coalescence of the 
hilltops is shown clearly in Fig.(^a). Why do those hilltops coalesce and how do they move? 
Will hilltop B in Fig.(P) move towards point A or point C? Answers to these questions can 
be addressed with the equation of motion of the hilltops. 

The curvature of the probe at time t is equal to l/(Xt). The probe becomes broader 
and broader with time. As long as the absolute value of the local curvature around the 
contact points, A or C, is much larger than the absolute value of the curvature of the 
probe, the positions of these two contact points do not change. In this case, the height 
difference h(xA,t) — h(xc,t), where xa and xc are the horizontal coordinates of points A 
and C respectively, is also a constant. It can be written as 

h(x A ,t) - h(x c ,t) = 

(h(x A , t) - h(x B , t)) - (h(x c , t) - h(x B , t)). 

The curve between A and B and the curve between B and C of the new surface are two 
parabolas. With h(x\,t) — h(xB,t) = (x A — £b) 2 /(A£) and h(xc,t) — h(xB,t) = (xc — 
XB) 2 /(Xt), we obtain the equation of motion of the hilltop, 

dx B ! ( x A + x c 



This equation implies that the hilltops always move toward the closest of the two contact 
points. Moreover, the closest contact point is also the higher one. If two neighboring hilltops 
share the same contact point, and that contact point is the highest one for both hilltops, 
then these two hilltops will merge into one around the position of the shared contact point 
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as illustrated in Fig.(|^). A sequence of coalescence events of the hilltops forms a tree-like 
structure of world lines as in Fig (|7|a). 

The dynamics of the coalescence of the hilltops in the noiseless KPZ equation is totally 
deterministic and only depends on initial surfaces. In the noisy case, the dynamics be- 
comes stochastic, however, the qualitative behavior of the coalescence of the hilltops is still 
sustained even when the surface is perturbed with randomly depositions. The nonlinear 
term causes the coalescence of the hill tops and remains essential even when the dynamics 
is stochastic. To check that the nonlinear term is responsible for the coalescence in the 
stochastic dynamics, we apply our upwards moving PRW model to EW-type surfaces. The 
upwards moving PRWs for EW type surface do not coalesce although the density of the 
PRWs is slightly higher on the hilltops. 

VI. DOWNWARDS MOVING PASSIVE RANDOM WALKERS 

So far we have focused only on the case in which k and A have opposite signs, i.e., when 
the PRWs move upward(/t > 0) in the KK model (A < 0). The discussion in the preceding 
section suggests that the shape of the surface has very different characteristics between a 
hilltop and a valley bottom on the surface with noiseless KPZ dynamics. Namely, hilltops 
are sharp with discontinued slopes while valley bottoms are rounded and with continuous 
slopes. The symmetry is broken by the nonlinear term. We also found that valley bottoms 
vanish by themselves, instead of coalescing with others valleys. 

If k is negative, the PRWs move toward the valley bottoms. It is interesting to see how 
the PRWs respond to such qualitative aspects of the surface in the presence of stochastic 
noise. In the deterministic case, the hilltops and valley bottoms intertwine, the number of 
valley bottoms and the number of hilltops both decrease. This is not true for the stochastic 
case. New microscopic hilltops and valley bottoms are being formed constantly by the noisy 
depositions of particles. Is the coalescence of the PRWs stable under such perturbations? Of 
upwards moving PRWs in the KK model, coalescences are stable against the noise. As we 
will see next, the coalescence of the downwards moving PRWs is not stable against random 
depositions. 

Consider the upwards moving PRWs. Suppose noise splits the hilltops by creating a small 
valley between them. This small valley is unstable against the dynamics as we have seen in 
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the preceding section. These two hilltops will merge again soon as a result of the nonlinear 
term in the KPZ equation. The nonlinear KPZ dynamics stabilizes the coalescence of the 
upwards moving PRWs. 

Next, imagine several downward moving PRWs moving into the same valley. The aggre- 
gation of the PRWs is not stable against the noisy perturbation. Depositions of particles 
on the valley split the valley. A PRW on the original valley bottom can go to either one 
of those two valleys with equal probability, so both valleys will be occupied by PRWs. Do 
the new created subvalleys always merge back into each other? A valley bottom does not 
move in the noiseless case; they will remain separated by the hilltop. The hilltop does not 
vanish unless it is driven away to merge with other hilltops, which is a much less likely 
process than the splitting. Therefore, dynamics like this does not stabilize the aggregation 
of downwards moving PRWs. They tend to be separated under the growth dynamics and 
we expect that the dynamics of downwards moving PRWs, trapped in valley bottoms, is 
dominated by the perturbation of random deposition rather than the deterministic nonlin- 
ear term in the dynamics. In Fig.([|c), we show a typical simulation for downwards moving 
PRWs. Instead of coalescence, we find that the aggregation clusters of PRWs are not stable 
and split constantly. 

Coalescence of the upwards moving PRWs for negative A is an important feature of 
KPZ-type surfaces. Not only does it ensure that the PRWs end up trapped on the global 
maximum; it also reveals the space-time structure of KPZ-type surfaces in detail, i.e., hilltops 
(valley bottoms) coalesce as time evolving for negative (positive) A. 

Our analysis on the KPZ equation applies to the phenomena of coalescence of the shock 
waves in the inviscid(z/ — > 0) noiseless Burgers equation. The velocity in the Burgers equation 
is the slope on the KPZ-type surfaces, i.e., v = d x h. By differentiation of eq. (|l|), it becomes 

dv . dv d 2 v . . 

¥ = Al "& + &3 + " B ' (20) 

The noise term satisfies (t]b{x, £)t7b(0, 0)) ~ d%.5(x)6(t). Because the slopes of the surface 
are discontinuous at the hilltops, the hilltops on KPZ-type surfaces correspond to the shock 
waves in the Burgers equation. The coalescence of the hilltops of the surface indicates 
that the shock waves coalesce, as illustrated in Fig. (0(b)). Coalescences of shock waves 
described by the noiseless Burgers equation have been reported in a numerical study of 
inelastic collisions of particles in 1+1 dimensions Bohr and Pikovsky also showed that 
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the zeros in the velocity field of the Burgers equation coalesce |TT|. The above analysis 
for the coalescence of KPZ hilltops translates directly into the coalescence of shock waves 
in the Burgers equation with noise. Our analysis on coalescences of PRWs on a growing 
surface presented here not only provides further intuitive and detailed understanding of the 
mechanism of these coalescence phenomena but also explains, in both deterministic and 
stochastic dynamics, the different behaviors between hilltops and valleys, which are distinct 
types of zeros in the velocity field in the Burgers turbulences. 

VII. RIVER-LIKE NETWORK AND PASSIVE SCALARS 

To further understand the stability of the aggregations of PRWs, we study how random 
forces, in addition to the driving forces from surface slopes, are applied to the PRWs directly 
and change the tree-like structure of the PRW world lines in Fig.(^a). Is the tree-like 
structure stable? Random noise modifies the equation of motion of the PRW as 

Qu d 

= K—h(x,t)\ x=u{t) +r] w (x,t), (21) 

with r] w the uncorrelated noise, (r) w (x, £)t/ w (0, 0)) = D w 5(x)5(t). D w is the diffusion constant 
of the PRW if k = 0. 

First consider a single PRW on a surface. The second term on the right hand side 
of eq. (|2lD introduces the ordinary diffusive behavior, i.e., Au(t) ~ (D w t) l l 2 . The scaling 
behavior due to the first term, Au(t) ~ (i£ w £) 1 / Zw , (K w is a constant determined by n and 

A) still dominates the large scale behavior since l/z w > 2. However, the PRW appears as 

i 

diffusive instead of super-diffusive at length scales smaller than l c ~ (^ SL ) 2_Zw • 

The diffusive noise also affects the coalescence of multiple PRWs. For small perturbations, 
the aggregation is stable because the PRWs cannot escape from the hilltop by diffusion if the 
hilltop is high and large enough. The net effect of diffusion simply broadens the coalescence 
cluster. However, if we increase the strength of the random noise, the PRWs can escape 
from the hilltop and be caught by other hilltops nearby with finite probability. In this case, 
the cluster of coalescence PRWs splits and the world lines of the PRWs form a braid-like 
network instead of a tree structure. A typical configuration of such a braid-like network is 
shown in Fig.(||b). Both the tree-like and the braid-like network resemble river-like networks 



of different length scales in Nature[16 
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The crossover scale between tree-like networks and braid-like networks is given by l c . We 
observe diffusive behaviors locally at scales less than l c . The braid-like networks emerge at 
the crossover scale, l c , and the tree-like structures are recovered at scales much larger than 

Furthermore, we can consider the evolution of the probability distribution function, 
P(x, t), of finding a particle at position x at time t. Note that the total number of PRWs is 
conserved in our model. Therefore, P(x, t) has the standard conserved form, 

^ = V.;(:M), (22) 

where j(x,t) is the PRW current. The current is the sum of the advected part and the 
diffusive contribution, i.e., j(x,t) = KvP(x,t) + D w d x P(x,t). Hence, we have 

= K V(vP(x, t)) + D w V 2 P(x, t). (23) 

at 

This is simply the equation of motion of passive scalars, P(x,t), in a fluid. In the limit 
D w — > 0, the passive scalars will be concentrated at a globe hilltop as we have seen in the 
previous sections if k and A have opposite sign for a finite system. 



VIII. CONCLUDING REMARKS 



Motivated by the domain wall dynamics on nonequilibrium fluctuating surfaces, we study 
the dynamics of passive random walkers on KPZ-type growing surfaces. We show that, 
although the coupling rule between a PRW and a surface is defined locally, the PRW typically 
reaches the maximum of the surface over a distance i in a time scale t ~ x Zs and "feels" 
the fluctuations of the surface over the same length scale. Two PRWs separated by Ax 
coalesce after At ~ (Ax) Zs . The fluctuations of the positions of such PRWs follow the same 
dynamical scaling as the KPZ fluctuations. We verify this scaling behavior numerically. 
Tracing the paths of the PRWs on KPZ-type surfaces is an effective method to measure the 
dynamical exponent directly in the stationary state. 

In addition to the scaling, the dynamics of passive random walkers reveals detailed specific 
space-time structures of a KPZ-type growing surface, i.e., hilltops coalesce with each other 
and the world lines of the PRWs form self-affine tree-like structures. We provide an analytical 
argument explaining this phenomenon based on the noiseless Burgers-KPZ equation. We 
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FIG. 8: World lines of passive walkers for different parameters. The lines are the paths of 128 
PRW on a system of size L=128. The darkness is proportional to the number of PRWs for each 
site. In (a), the system is simulated in the limit l c = (-^p) 2_z —> for upward moving PRWs. 
In (b), we choose l c ~ L, the system size, (c) is the result of downwards moving PRWs in the 
limit l c —* 0. Even without a diffusion term in eq. (|2l]), the downwards moving PRWs do not 
coalesce. The dynamics of the downwards moving PRWs on KK-type surface is dominated by the 
uncorrelated random depositions. 
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show how the nonlinear term (d x h) 2 is responsible for this nontrivial coalescence phenomenon 
in both the noiseless and noisy case. The noiseless KPZ equation also allows us to understand 
why downwards moving PRWs do not coalesce, due to the asymmetry between hilltops and 
valley bottoms. 

We want to emphasize that the effectively attractive interactions which makes the PRWs 
coalesce are not only topological but also robust against perturbations. For KPZ-type 
growth, the (d x h) 2 term breaks the symmetry between hilltops and valley bottoms during 
growth and introduces nontrivial dynamics of the hilltops for negative A and the valley 
bottoms for positive A. The singularities on KPZ-type surfaces generated by the nonlinear 
terms stabilize the aggregations of the PRWs against other perturbations. The coalescence 
of PRWs reveals one of the important features of the complicated space-time structure of 
the surface. An interesting extension of this study is to see how to generalize this to different 
universality classes of surface dynamics. We already point out that in EW growth, because 
of the particle-hole symmetry, PRWs do not coalesce. For other nonlinear models, in which 
the nonlinear terms break up-down symmetry, PRWs may be driven in other nontrivial ways 
and the surface may show interesting space-time structures. 

Depending on the strength of the noise applied to the PRWs, the world lines of PRWs 
form tree-like or braid-like networks in space-time. The tree-like or the braid-like networks 
resemble river type networks in nature. The braid-like networks are the result of two com- 
peting mechanisms, namely, the coalescence between PRWs caused by the fluctuations of 
surfaces and the diffusion of PRWs caused by directly applied noise. These two mechanisms 
define a crossover length scale Z c , and give rise to the braid- like networks at this scale. One 
might expect to see similar crossover phenomena in the study of passive scalars in Burgers 
turbulence. 

Previous studies on fluctuating nonequilibrium surfaces were focused on the scaling be- 
havior. This study shows that, beyond global scaling laws, such as the scaling of the global 
interface width, nonlinear surface growth dynamics leads to intriguing detailed aspects of 
their space-time structures. We hope studying these different perspectives of nonequilibrium 
dynamics can lead to a better understanding of these complex systems. 
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